Driven superconducting proximity effect in interacting quantum dots 
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We present a theory of non-equilibrium superconducting proximity effect in an interacting quan- 
tum dot induced by a time-dependent tunnel coupling between dot and a superconducting lead. The 
proximity effect, that is established when the driving frequency fulfills a gate- voltage-dependent res- 
onance condition, can be probed through the tunneling current into a weakly-coupled normal lead. 
Furthermore, we propose to generate and manipulate coherent superpositions of quantum-dot states 
with electron numbers differing by two, by applying pulsed oscillatory variations to the couplings 
between dot and superconductors. 
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I. INTRODUCTION 

The possibility to induce superconducting correla- 
tions in a quantum dot by the proximity of super- 
conducting leads has received substantial interest in 
the last decade. 1 ' 2 A variety of transport phenom- 
ena in hybrid normal/superconducting/quantum-dot sys- 
tems have been addressed experimentally, including the 
Josephson effect, 3 7r-junction behavior, 4 ' 5 superconduct- 
ing quantum interference, 6 Cooper pair splitting, 7 9 
and Andreev-level spectroscopy. 10 ' 11 Theoretical stud- 
ies have focused on the interplay of superconductiv- 
ity with the Kondo effect, 12 17 nonequilibrium Andreev 
transport, 18 25 and multiple Andreev reflection. 26 28 

In equilibrium, strong Coulomb repulsion in a quantum 
dot suppresses the formation of a pair amplitude, while 
different-time superconducting correlations are induced 
by higher-order tunneling. As a consequence, the Joseph- 
son current through a quantum dot tunnel-coupled to two 
superconductors will be at least of second order in the 
tunneling rate r. 29 ' 30 Experiments with quantum dots re- 
alized in carbon nanotubes 3 ' 5 ' 6 and InAs nanowires 4 are 
consistent with this scenario. However, it has been pro- 
posed that a non-equilibrium occupation of the quantum 
dot's states, induced by a voltage-biased normal lead, 
can lead to a finite pair amplitude on the dot. 31 ' 32 In 
this Article, we propose an alternative way, namely us- 
ing tunnel couplings that oscillate in time, to establish a 
non-equilibrium situation that supports superconducting 
correlations in the dot. 

Transport through few-electron quantum dots can of- 
ten be described by taking into account only a single 
orbital level. In this case, a finite pair amplitude in the 
dot describes a coherent superposition of this level being 
empty, denoted by |0), and doubly occupied \d). It can 
be generated by driving the system via oscillatory tunnel 
couplings. The resonance condition under which the pair 
amplitude becomes maximal involves the oscillation fre- 
quency and the energy difference S between the \d) and 
|0) states. The superconducting proximity effect can be 
probed by measuring the current into an additional, nor- 



mal lead weakly coupled to the dot. Finally, we show 
that suitable pulsed oscillatory variations of the system 
parameters can produce any coherent superpositions of 
the |0) and \d) states. 



II. MODEL 

We consider a single-level quantum dot tunnel cou- 
pled to superconducting and normal leads, which can 
be described by the Anderson-impurity model with the 
Hamiltonian H = Hqd + J2 v =s,n(^v + ^,77) • We de- 
note with d^ a (d a ) and c^ ka (c^ka) the creation (annihi- 
lation) operators for electrons with spin a =t,i in the 
dot and for single-particle states with quantum number 
k in lead 77, respectively. The quantum-dot Hamiltonian 
reads Hqd = ^2 a edld a + Ud^d^d^d^, where e and U 
are the single-particle energy and the Coulomb energy 
for double occupation, respectively. The leads are de- 
scribed by the BCS Hamiltonian H v = ^2 kcr ^k^ka^ko- — 
Ysk(\ c lkt c l)-kl + H.c.), witn the superconducting pair 
potential A^ = |A|e^ being non zero only for the su- 
perconductors. The coupling between dot and lead 77 is 
captured by the standard tunneling Hamiltonian Ht^ = 

Vv J2ka( c l]kcr^ + H.c). We assume the tunnel matrix 
elements to be independent of k and <r, and define the 
tunnel-coupling strength as = 27r\Vrj\ 2 ^2 k S(uj — e&). 

Since we focus on the proximity effect, whose micro- 
scopic origin is sub-gap Andreev reflection, we consider 
the limit of large superconducting gap, |A| — >> 00. In this 
limit, the tunnel coupling to the superconducting leads 
can be taken into account by an effective dot Hamilto- 
nian # QD ,eff = H QD - A e ff(t)d|d| - A* ff (t)d^d t , that 
shows a BCS-like term A e ff, whose strength is deter- 
mined by the tunnel couplings to the superconducting 
leads. The effective-Hamiltonian formalism has been de- 
rived and employed in the time-independent case, 32 ' 33 
but it remains valid in the time-dependent case since 
non-instantaneous contributions play a role only for fre- 
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FIG. 1. (Color online) Schematic setups. The oscillation in 
the effective coupling between dot and superconductor can be 
realized by (a) a time-dependent phase difference between the 
superconductors or (b) varying the tunnel-coupling strength 
between dot and superconductor. 



quencies larger than the superconducting gap |A|. 

We propose two different schemes to generate a time- 
dependent coupling A e ff(t). In the first scheme, shown 
in Fig. 1 (a), the dot is attached to one normal lead 
and two superconducting leads with a time-dependent 
phase difference = 2Qt. The latter can be created by 
monotonically increasing the magnetic flux enclosed by 
the superconducting loop or by directly applying a volt- 
age V = HQ/e between two superconductors. Assuming 
a symmetric capacitive coupling between the dot and the 
two superconductors, the voltage drop and, subsequently, 
the time-dependent phases of the superconductors are di- 
vided symmetrically, (/>i,2 = such that 



A eff (t) = -(r 5 ie int + r 5 2e 



-iQt\ 



(i) 



The two tunneling rates Tsi,2 can, in general, differ from 
each other. In the second scheme, shown in Fig. 1(b), 
the quantum dot is coupled to one normal and one super- 
conducting lead with Ts(i) = f s + ^I^ cos(Qt), yielding 
A e ff(t) = Ts(t)/2. In the following, we formulate the re- 
sults for the first scheme. The second scheme is covered 
by substituting Tsi and Ts2 with 5Ts/2. 

In order to obtain an effective description of the dot, we 
trace out the normal lead's degrees of freedom. This leads 
to the reduced density matrix of the dot pqd with matrix 
elements Pa(t) = (a\pQv(t)\f3) with respect to the basis 

states |0),| t)> I \)i \d) = ^|^||0) for the empty, singly- 
occupied with spin up or down, and doubly-occupied dot. 
The corresponding energies are Eq = 0, E± = = e, 
and Ed = 2e + /7 = 5, respectively. Among the off- 
diagonal elements, only the dynamics of Pq and are 
coupled to the occupation probabilities. The dot's pair 
amplitude is described by Pq (t) = (di(t)d^(t)). The time 
evolution of the reduced density matrix of the dot after 
time 0, at which the driving is switched on, is governed 



by the kinetic equation, 
d 



<ft i^(t)+t[ff Q D, eff (t),p QD (t)]| 



f dt i w${t,t i )p$;{t i ), (2) 
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where Wpj$, (t,t f ) are transition rates due to the coupling 
to the normal lead. These rates can be evaluated by 
means of a real-time diagrammatic technique. 32 

The evaluation of the diagrams is considerably sim- 
plified by employing the rotating- wave approximation, 
which is justified in the limit \Q\ ^> |A e ff(t)|. It con- 
sists of neglecting the non-resonant component, i.e., of 
replacing A eff (t) -+ \T S ie mt or A eff (t) -> \Ys2e~ int 
if 5 and have the same or opposite sign, respectively. 
In the following, we write all formulas for the first case. 
The second one is obtained by replacing Tsi with Ts2 
and Q with — O. Within in the rotating- wave approxi- 
mation the time evolution due to i^QD,eff, which enters 
the calculations of the rates Wqq, (t,t f ), takes the form 

U(t,t>) = e- i < t -*')(^|a)( C 7| +a+|0>(0| +a.\d)(d\ 



/3 + |0)(d|-/3-M>(0| 



7=± 



5-tt 



3 ±i£ 7)=F (t—t f ) 



7=± 



(3) 
(4) 



(5) 



Here £~ 



{ill 

energies with 7, 7' = ±1 and 2ca = y/(S — Q) 2 

By explicit calculation, we find that the kinetic equa- 
tions for P|, and P^+Po decouple from those for Pq = 
(P^) andP^ — Pq. The latter are conveniently expressed 
in terms of an isospin I = (RePg , IuiPq , (P^ — Po)/2). To 
first order in T^v, which describes the regime of weak tun- 
nel coupling between dot and normal lead (Tn <C feT), 
the isospin dynamics is governed by 



Q)/2 are Andreev addition 



r 2 
1 si- 



dt 



I(t)-I(t)xB<°»(t) = 



T N I dt' 
lo 



I(t') x B - R I(t') + A 



(6) 



Here, B^°\t) = (T S i cos(ftt), T S i sin(Ot), -S) is the ef- 
fective magnetic field inducing a fully-coherent evolution 
of the isospin in the absence of the normal lead. A finite 
tunnel coupling to a normal lead (r.h.s. of Eq. (6)) intro- 
duces a correction to the isospin rotation about vector 
B(t, t f ), a relaxation matrix R(t, t r ) and a vector A(t,t') 
associated with the generation of the isospin. Their (real) 
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components are determined by, 



B X A 


-iBy = 




(7) 




+ iB z -- 


= j ^[i-K-aD/MK^'), 
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= Bt J 


r ^[l + (a + -a*_)/H]e-^-*'), 
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(9) 


A x H 


-iAy = 


J ^( / 3 + + ^)/( w)e Mt-0 ) 


(10) 


A z = 


-Bef 


^[(a + + al)/H-l]e-^-*'). 


(11) 



Here f(oj) = 1/[1 + exp(cj//c#T)] indicates the Fermi 
function. We note that the relaxation matrix is diago- 
nal and R yy = R xx 

In the absence of a normal lead, T^v = 0, Eq. (6) 
describes a fully coherent evolution of the isospin un- 
der B(°)(£). The solution in this case can be written 
as I(t) = 0(£)I(O) with Q(t) being a rotation matrix 
with time-dependent rotation axis and angle. To solve 
Eq. (6) in the presence of the normal lead we, first, 
define the auxiliary quantity J(t) — _1 (£)I(£), which 
describes the isospin in the time-dependent coordinate 
system that follows the coherent dynamics induced by 
B(°)(t). The kinetic equation for J(t) is still an integro- 
differential equation. However, when restricting to first- 
order tunneling processes, we can substitute J(t) for J(t') 
in the integrand. 34 In fact, the master equation shows 
that the derivatives appearing in the general expansion 
J (0 = E£Lo(* - t') n [d n J(t)/dt n ]/n\ are at least of first 
order in Tj^. Plugging them into the integrand on the 
r.h.s. of the master equation for J(t), these terms corre- 
spond to at least second order in Tat and are, therefore, 
neglected. Another simplification we make is the replace- 
ment of the lower limit of the integral by — oo which is 
justified when t is much larger than the characteristic 
decay time of the integral kernel. Explicitly performing 
the integral leads to a pure differential equation for J(t) 
that can be solved analytically. This yields lengthy ex- 
pressions that we do not show here. The time evolution 
of I(t) is, then, obtained by going back to the original 
coordinate system. 



III. RESULTS 

In the following, we consider two different regimes: 
(1) ^> Tn and (2) Ts ~ IV- In the first case, when 
the dot is strongly coupled to the superconductor (s) 
and weakly coupled to the normal lead, the system 
is mostly governed with Rabi-resonance physics. It is 
in this regime, that coherent manipulations of the dot 
charge states, detectable by coherent oscillations in the 
current, can be achieved. In the second regime, the Rabi 
resonance still plays an important role, but features 
related to quantum stochastic resonance become relevant 
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FIG. 2. (Color online) Pair amplitude strength Aqd as a func- 
tion of oscillation frequency Q for different detunings 5. The 
parameters are T^i = lksT, Ts2 = and U = lOksT. 

Strong pairing is seen around the resonance points Q ~ ±6. 
The apparent discontinuity at Q = is due to the failure of 
the rotating- wave approximation at small frequencies. 

since the driving and relaxation have the same order. 



A. Strong coupling to the superconductor 

Here we investigate the regime in which the coupling 
of the dot to the superconductors is much stronger than 
to the normal lead (r^ ^> Tn). First, we discuss the 
steady-state behavior when the system reaches a constant 
J = Joo for times t^> 1/Yjsf. The corresponding isospin, 
loo(^) = @(^)Joo7 rotates about the z-axis with frequency 

i.e., in addition to the constant ^-component, there is 
an oscillatory in-plane part of the isospin. 
In the limit of large when /(£ 77 ') for all Andreev 
addition energies are approximately zero, the isospin can 
be written as 

The in-plane components of the isospin indicate super- 
conducting correlations, described by the pair amplitude 
(d^df) = Po = Ix + Hy = e* m A QD , where A QD charac- 
terizes the (time-independent) strength of the proximity 
effect. For large detuning, \S\ ^> Tsi, the z-component of 
B(°) tries to keep the isospin aligned with the z-direction, 
which corresponds to a negligibly small Aqd- However, 
once the driving frequency fulfills the resonance condi- 
tion \5 — 0| < Tsi, the in-plane component of the isospin 
becomes significant. If S and Q have different sign then 
the resonance condition is replaced by \S + fl\ < Ts2- 
Figure 2 shows Aqd as a function of driving frequency 
ft for various values of detuning S. Close to the reso- 
nance conditions and more precisely at Q — S = ±Tsi and 
Q + S = tFs2, the pair amplitude attains the maximum 
and minimum values of ±1/4 as one can see from the 
isospin solution given by Eq. (12). Exactly at resonance, 
Aqd displays an abrupt sign change, which is reminis- 
cent of the — 7r transition in the Josephson current. 4,32 
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FIG. 3. (Color online) Normal-lead current as a function of 
frequency Q for different values of 5. The parameters are 
Tsi lk B T, T S 2 3k B T, and U 10k B T. The current 
shows resonance peaks at f2 = +5 and —5 with widths Tsi 
and Ts2, respectively. 



Furthermore, we remark that for \S\ < U the pair am- 
plitude is strongly suppressed since in this regime single 
occupation of the dot is favored. 

The nonequilibrium proximity effect can be probed 
by measuring the current into a normal lead, tunnel- 
coupled to the dot. The current in lead N can be 
computed by means of the general formula 7/v(£) = 

Eaa'/3' Jo* dt,wN o$?' (*> 0^0'' (0> where the current rates 
W^j^ (£, t') take into account the number of charges that 
flow in or out of the normal lead. In our system, the cur- 
rent up to first order in Tn depends only on the isospin 
degree of freedom and is given by 



I N (t) = -T N 



f dt' 
Jo 



{R z 



A z + (IxB) z }. (13) 



In the steady-state limit the current turns out to be con- 
stant and for large frequency and detuning, using the 
isospin given by Eq. (12) we obtain the current with a 
Lorentzian form, 



I 



N 



ft 



N 



F 2 
1 SI 



+ -0)2 



(14) 



In Fig. 3 the current is plotted as a function of frequency 
for various values of 6 > 0. Resonance peaks with width 
T^i and Ts2 appear at ft = +5 and —5, respectively. For 
negative values of 5, the current flows in the opposite 
direction, and the position of the peaks with width Tsi 
and Ts2 are interchanged. 

In order to discuss the transient behavior of the driven 
dot we assume a positive large S > U which in equilib- 
rium, causes the dot to be empty, I = (0,0, —1/2). After 
switching on the resonant driving (Q = S) at t = 0, I z de- 
cays on the time scale 1/Tjv, accompanied by oscillations 
with frequency Tsi- Pulses of finite length in time allow 
to generate and manipulate coherent superposition of the 
empty and doubly-occupied states. For instance a pulse 
with duration t = n/(2Tsi) causes a rotation of isospin 
from z into the x-y-plane. In order to probe these Rabi- 
like oscillations via the (intrinsically stochastic) tunnel 
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FIG. 4. (Color online) The time-averaged current in the 
normal lead as a function of pulse duration t on . The inset 
shows the real-time variations of the dot charge when the 
driving is switched on and off periodically with durations 
ton = OAStt/Tn and t Q ff = 5/Tjv, respectively. The val- 
ues of the parameters used for this figure are U = 10k bT, 
Q = 6 = 50k B T, Tsi T S2 2k B T, and F N 0.2k B T. 



current into the normal lead, we propose to periodically 
switch the driving on and off with durations t on and t ff 
(where t Q R ^> I/Tat). The time-averaged current into 
the normal lead then shows damped oscillations as func- 
tion of t on and reaches the steady state value er^r/ft for 
ton ^> t ff, see Fig. 4. The underlying real-time dynamics 
of the dot charge Q = — e(l + 21 z ) shows damped Rabi- 
like oscillations (pure decay) when the driving is switched 
on (off) as illustrated in the inset of Fig. 4. 



B. Weak coupling to the superconductor 

We now turn to the regime of weak coupling to both 
superconducting and normal leads (r^ ~ Tn <C fc^T). 
For a systematic perturbative treatment of the isospin 
dynamics, we count both Ys and Tn but also the de- 
tuning from the resonance condition 5 — Q as one or- 
der in the perturbation expansion. As a consequence, 
we drop all higher-order terms of Ts and S — Q in the 
master equation. This can be easily achieved by the re- 
placement £ 11 >( k €a 0) = (j'U + Q)/2 for the Andreev 
addition energies. As a result the form of the kinetic 
equation (6) for the isospin remains unchanged but B, 
R, and A are modified according to the replacement of 
Andreev addition energies. In addition the effective mag- 
netic field B(°) (zeroth order with respect to T^) should 
now be considered as sum of a zeroth-order, (0,0, —J), 
and a first-order part, (Tsi cos(Ot), Tsi sin(Qt), 0), with 
respect to Ts- Therefore in this case the matrix 6 used 
in the definition of auxiliary quantity J = corre- 
sponds to the coherent evolution under (0,0, —5). 

Despite these modifications, the method to obtain the 
isospin dynamics is the same as in the strong coupling 
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FIG. 5. (Color online) Real and imaginary components of 
the pairing strength Aqd as a function of frequency SI for 
various values of tunnel couplings Tn and Tsi> The inset 
of the lower panel shows the dependence of the imaginary 
part of Aqd on the relative value of the tunnel couplings 
IV/Isi assuming SI = 5 and Tsi = O.lfcsT. The values of 
parameters used in this figure are U = 10k bT and 5 = 2U. 
Here the pairing strength around SI ~ 6 is shown and for 
negative values SI ~ — S the plots can be obtained by mirror 
reflections with respect to and with widths determined by 
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case. For the steady- state limit, we find 
A/K 
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in which ip is the digamma function. Here, we follow the 
same convention for the rotating-wave approximation as 
before: only the results for the case of 5 and SI having the 
same sign are presented. The results for 5 and SI having 
opposite signs can be obtained by the replacement SI — » 
— SI and T^i — » Ts2- From the isospin form given above, 
we find an oscillatory pair amplitude with frequency SI 
and a complex- valued strength 
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r 2 

1 SI 



(5-Sl-B) 2 ' 



(19) 



Similar to the strong-coupling limit, the pairing strength 
has a resonant form as function of SI. There are, how- 



ever, differences. First, Aqd now depends on both cou- 
plings Tsi and Tn and, second the resonance point is 
slightly displaced by an effective exchange field B. As 
one can see from Fig. 5 the real part of the pairing 
strength behaves in a similar way to Aqd in the strong- 
coupling limit and passes through a maximum/minimum 
at S — SI — B = ±y/lZ 2 + T 2 S1 accompanied with a sign 
change at 5 — SI — B = 0. On the other hand, the 
imaginary part, which was absent in the strong-coupling 
regime, shows a single resonance peak at S — SI — B = 0. 
Interestingly, although the maximum value of the real 
part monotonically decreases by increasing IV, the peak 
value of imaginary part has a maximum for = ^si- 
This is clearly seen in the inset of Fig. 5 and can be 
understood as a manifestation of the so-called quantum 
stochastic resonance, 35,36 a phenomenon which appears 
when the response of the system to a driving signal is 
increased by a certain level of dissipation 37 or noise. 38 
We note that away from the resonance point, the above- 
mentioned signature of the quantum stochastic resonance 
disappears: the absolute value of the pairing, which is 
physically more relevant, does not need an optimal value 
of Tn to become large. As one can easily check from 
Eq. (19), | Aqd | passes through a maximum \ A\/21Z when 
(5 — SI — B) 2 +7Z 2 = T 2 S1 . This is not achievable for large 
dissipation Tn ^> Tsi, but always possible for T si > IV, 
with an optimal pairing amplitude of 1/4, in accordance 
with the result for the strong-coupling regime. 

To obtain the current we use (13) and introduce the 
same modifications as the for kinetic equation. We plug 
in the isospin solution (15) and get the steady-state cur- 
rent, 



in 



2T 2 S1 A 



nn 2 



■rli 



(5 -SI- B) 2 



(20) 



which, again, has a Lorentzia n form wit h displaced reso- 
nance S - SI = B and width y/K 2 + T 2 S1 . For Tsi > IV, 
the current reaches the value Eq (14) of the strong- 
coupling regime. 

We end up this subsection by commenting on the tran- 
sient behavior in this regime. Since we are assuming weak 
coupling to the superconductor T$ ~ IV, the relaxation 
mechanism dominates over the coherent evolution pro- 
duced by the driving and no coherent oscillations can be 
seen. This means that a coherent-manipulation scheme 
of the charge states requires a stronger coupling to the 
superconductors than to the normal lead. 

Finally, we comment on possible experimental verifi- 
cation of the predicted effects. Considering typical ex- 
perimental values ksT ~ 1 /ieV, |A| ~ 0.15 meV (for 
Al/Ti) 4 , S and SI must be on the order of 10 jaeV to fulfill 
the requirement T N <C k B T,T s <C \S\, \Sl\ <C |A|. There- 
fore, for the setup where magnetic flux is used to pro- 
duce a time-depedent tunnel coupling, a magnetic field 
variation rate dB/dt ~ 10 5 T/s is necessary for a super- 
conducting loop with area 100 /im 2 . 
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IV. CONCLUSIONS 

We have proposed the possibility to induce supercon- 
ducting correlations in a single-level quantum dot via 
a time-dependent tunnel coupling to a superconduct- 
ing lead. A finite pair amplitude can be generated if 
a resonance condition for the driving frequency in rela- 
tion to the energy difference between empty and doubly- 
occupied dot is fulfilled. The existence of this nonequilib- 
rium proximity effect can be probed by measuring a finite 
steady-state current into a weakly-coupled normal lead. 
Furthermore, Rabi-like oscillations between empty and 



doubly-occupied dot states may be probed in the tunnel 
current by applying periodically pulsed oscillations of the 
tunnel coupling to the superconductor with tunable pulse 
duration. 
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